Decoherence of many— body systems due to many— body interactions 
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We study a spin-gas model, where Ns system qubits are interacting with Nb bath qubits via 
many-body interactions. We consider multipartite Ising interactions and show how the effect of de- 
coherence depends on the specific coupling between the system and its environment. For instance, 
we analyze the influence of decohenerce induced by fc-body interactions for different values of k. 
Moreover, we study how the effect of decoherence depends on the correlation between baths that are 
coupled to different individual system qubit and compare Markovian with non-Markovian interac- 
tions. As examples we consider specific quantum many-body states and investigate their evolution 
under several different decoherence models. As a complementary investigation we study how the 
coupling to the environment can be employed to generate a desired multipartite state. 



I. INTRODUCTION 

In any realistic secnario the system of interest is never 
completely decoupled from its environment. This cou- 
pling causes the system to decohere vifhich is believed to 
be the main reason why macroscopic quantum systenrs 
do not occur in nature. The effect of decoherence has 
widely been investigated in the literature. There, dif- 
ferent decoherence models have been assumed based on 
interactions with a harmonic oscillator bath [1], a spin 
bath [2] or a spin gas [3] . 

We will consider here a microscopic decoherence 
model, where both, the system and the environment are 
described by qubits. For a general spin model, where Ns 
system (S) qubits are coupled to Nb bath (B) qubits, the 
interaction would be described by a general unitary, Usb, 
acting on the whole system. Assuming that the initial 
state of the bath would be described by the state ps , the 
evolution of the system qubits would be governed by the 

completely positive map £{ps) = ^"csiUsBPs (X" PbUsb)- 
Due to the exponential scale of the dimension of the 
Hilbert space the investigation of the evolution of the 
system qubits under such a general map is unfeasible. 
Thus, one has to restrict the coupling Usb to a certain 
class of unitaries and therefore to a certain class of cou- 
pling Hamiltonians Hsb- A reasonable choice is to con- 
sider only generalized Ising interactions, which describe 
the effect of multi-qubit collisions. The coupling Hamil- 
tonian for this choice Hsb is of the form Hsb = X^i Q^i^'L 
with i denoting a {Ns + NB)~h\% string. Each term in 
the Hamiltonian corresponds to a certain collision. The 
unitary Usb ~ g-i-Hsst^ which governs the evolution of 
the state, is a phase gate, i.e. a unitary which is diag- 
onal in the computational basis. Due to the fact that 
all terms in the Hamiltonian commute, Usb can be writ- 
ten as a product of phase gates, where each factor cor- 
responds to a certain collision. The order of the phase 
gates, i.e. the number of qubits on which the phase gates 
acts non-trivially, is equivalent to the number of qubits 
which collided. The value of the phase in the phase gate 
depends on the interaction strength. In this article we 



will describe the evolution of the system in terms of phase 
gates for arbitrary phases. That is, we will not consider 
the evolution as a function of time, but will rather as- 
sume that due to a certain physical interaction a specific 
phase gate, Usb is given. The investigation of such a 
general model is only feasible because all the phase gates 
corresponding to different collisions commute. 

In this article we investigate first the effect of decoher- 
ence induced by multi-qubit Ising interactions. Second, 
we will show how the coupling between the environment 
and the system can be employed to prepare the systenr 
qubits in a desired state. The first part of this article 
constitutes a generalization of the work presented in [4] , 
where the decoherence of a spin gas has been investi- 
gated. In contrast to [4] we do not restrict ourself to 
bipartite interactions, but consider multipartite interac- 
tions and the resulting decoherence. Obviously, consid- 
ering many-qubit phase gates leads to a much bigger va- 
riety of different interactions. For two-qubit phase gates, 
the only relevant interaction is the one where one bath 
qubit is interacting with one system qubit. In the case of 
three qubit gates, we could for instance consider 3-qubit 
interactions between either 2 system qubits and 1 bath 
qubit or between 1 system qubit and 2 bath qubits. Apart 
from that, considering several system qubits, which are 
interacting for instance with the bath qubits via a three- 
qubit phase gate, the bath qubits could either be inde- 
pendent, or not. In this case, non, one, or both of the 
bath qubits could interact with both system qubits. The 
effect of decoherence will strongly depend on the corre- 
lations within the bath. Considering more general phase 
gates (acting non-trivially on k systems) allows for cor- 
respondingly more different kind of interactions. 

The aim of this paper is to get a better understanding 
of the effect of non-local decoherence on a multipartite 
system. This will be achieved by comparing the influ- 
ence of decoherence for different interaction models. For 
instance, we will compare the effect of two-qubit phase 
gates to the one of three and more qubit gates. Moreover, 
we will compare Markovian to non-Markovian processes. 
We will furthermore investigate how the correlation in 
the bath qubits affect the decoherence. Even though we 



will consider the most general dephasing maps, our main 
focus will be on those maps, where each system qubit is 
only interacting with bath qubits, i.e. there is no inter- 
action between two, or more system qubits and the bath. 
We will call these maps purely dephasing maps. The rea- 
son for focusing on purely dephasing maps is that they 
correspond to the realistic situation where the processes 
in which several system qubits collide with bath qubits 
can be ignored. Another reason for focusing on purely 
dephasing maps is that they cannot generate entangle- 
ment. Thus, they are perfectly suited to investigate the 
effect of decoherence on the entanglement and coherence 
of the system qubits. 

The outline of the paper is the following. First of all 
we will review the notion of a certain class of multipar- 
tite states, called locally maximally entangleable (LME) 
states. As we will see, the reduced states of LME states 
(LMESs) play an important role in our considerations 
of dephasing maps. The projected-entangled-pair state 
(PEPS) picture will be used to compute their reduced 
states. In Sec. Ill we will consider the most general de- 
phasing maps and will derive a simple relation between 
the maps and the corresponding LMES. After that, we 
will focus (apart from Sec. V) on purely dephasing maps 
which result from interactions where each system qubit 
is only coupled to baths qubits and not to other system 
qubits. We will show that all these maps are separable, 
i.e. they can be written as a convex combination of local 
maps and will express them in terms of simple local maps. 
In order to study how the bath-correlation influences the 
decoherence, we will derive a simple expression of the 
maps which describe the following scenario. We consider 
m bath qubits which are all interacting with each system 
qubit. Each of the system qubits is also interacting with 
some other bath qubits, which are coupled to the m bath 
qubits, but not to any other system qubit. The influence 
of correlated noise will be investigated by considering the 
maps described above for different values of m. In Sec. 
IV we will analyze how the effect of decoherence depends 
on the various coupling scenarios. In particular, we will 
investigate the Markovian case and compare the effect 
of decoherence to the non-Markovian case for a general 
input state (Sec. IV A). Furthermore, we will compare 
the effect of phase gates of order k for different values of 
k. As examples we will determine the effect of decoher- 
ence on arbitrary single system qubit state, two system 
qubits states and certain multipartite states (Sec. IV B). 
We will show that the intuition, that the entanglement 
shared by the system qubits is less disturbed, if the baths 
of the system qubits arc stronger correlated, that is m is 
large, is not true. It will be demonstrated that for cer- 
tain multipartite states the entanglement is more robust 
under the decoherence of independent baths than under 
correlated ones. In Sec. V we will consider the com- 
plementary process to the one considered above. There 
we will employ the interaction between the system and 
bath qubits in order to prepare the system qubits in a 
desired multipartite state. In particular, we will show 



how LMESs can be generated via unitary interactions 
between the system qubits and the bath qubits. 

Through out the paper we will use the following no- 
tation. The subscript of an operator will always denote 
the system it is acting on, or the system it is describing. 
We denote by i the classical bit-string (ii, . . . ,i„) with 
ik & {0, 1} Vfc g {1, . . . , n} and |i) = |ii, . . . , i„) denotes 
the computational basis. Normalization factors as well 
as the tensor product symbol will be omitted whenever 
it does not cause any confusion and 1 will denote the 
normalized identity operator. 



II. LOCALLY MAXIMALLY ENTANGLEABLE 

STATES 

Recently a new classification of multipartite states 
has been proposed in order to study the properties of 
pure quantum states describing n qubits [5]. There, 
the following question has been addressed. For which 
states, 1^), do there exist local controlled operations, 
^i — Si ^i ^ N)/ (*l I where [// are unitary operations 
acting on system I and |i)j (j| is acting on the auxil- 
iary system attached to system I, such that the state 
Ci «) C2 «) . . . «) C„ |*)g (|0) + |1))®" is maximally entan- 
gled between the system (S) and the auxiliary system 
(A). States for which this is possible were called Locally 
Maximally entangleable States (LMESs). Important ex- 
amples of these states are all stabilizer states. In [5] it 
has been shown that any LMES is (up to local unitary 
operations) of the form 



1^) = \ — 
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(1) 



where a-i € R and t/^^ denotes the diagonal unitary oper- 
ator with entries e'"' (phase gate). Note that the 2" real 
phases ai characterize the LMES (up to local unitaries) . 
From Eq (1) it can be easily seen that any LMES can be 
prepared by applying generalized phase gates, acting on 
up to all qubits, to a product state, i.e. non-entangled 
pure state. That is, any LMES |^) can be written as 



I*) 



Ul,...,nX{U,, 



Wu. 



(2) 



where C/^^ ,...,ik is a phase gate acting non-trivially on I 
qubits. For instance, L/123 = 1- (l-e"*!^^) [in^ (lll|i23> 
with 0123 € 1? is a three-qubit phase gate. It is straight- 
forward to see that in this hierarchical way the 2" phases 
ai can be generated. Thus, any LMES can be prepared 
using generalized phase gates, which could result from 
a generalized Ising interaction. If a(i) is a polynomial 
of degree k (as a function of i = (ii, . . . , i„)) then the 
corresponding state can be prepared using only A: -body 
interactions. Thus, the correlations in the coefficients are 
directly related to a preparation scheme and therefore to 
the entanglement contained in the state. 



Apart from those properties, LMESs have the fohowing 
properties [5]: (i) According to their definition, LMESs 
are characterized by the fact that their global information 
can be washed out by maximally entangling the system 
qubits to local auxiliary qubits. (ii) A state is LME iff 
it can be used to encode locally the maximum amount 
of n independent bits, (iii) For any LMES, |^), one can 
construct a complete set of commuting unitary obscrv- 
ables such that |^) is the unique eigenstate with eigen- 
value one for all these observables (the so-called gener- 
alized stabilizer, see also Sec. V). (iv) LMESs can have 
an exponentially large quantum Kolmogorov complexity 
[6]. That is, for certain LMESs it requires exponentially 
many classical bits to describe the state. 

Note that this classification of multipartite states is in- 
teresting for several quantum information theoretic tasks. 
First of all, (i) is equivalent to the fact that all the infor- 
mation contained in a LMES can be coherently " copied" 
into local auxiliary systems by local unitary operations. 
This seems to be a crucial property shared by those states 
which are useful for quantum information tasks, like one- 
way quantum computing, quantum error correction and 
quantum secret sharing. Second, due to the preparation 
scheme [Eq. (2)], LMES can be entangled in many dif- 
ferent, but hierarchical ways. Product states but also 
stabilizer states arc all LME. Stabilizer states or more 
generally, weighted graph states [3, 4, 7] arc those LMES 
which require only two-qubit phase gates for their prepa- 
ration. LMES, however, can be much more entangled 
than the stabilizer states. Moreover, there exist, in con- 
trast to stabilizer states, complex LMES [(iv)], which are 
necessarily highly entangled [6] . 

In the following subsection we will consider the reduced 
state of LMESs, which will be relevant for our consider- 
ations of dephasing maps. In Sec. V we will show how 
LMESs can be generated using either dissipation or uni- 
tary interactions. 



A. Reduced state of LMESs 

In Sec. Ill we will show that dephasing maps arc 
strongly related to the reduced states of LMESs. Due 
to this fact we will present here some methods for com- 
puting them. 

Let us consider an arbitrary n-qubit LMES, |^) = 



U\ 



\(gm 



^i e'"' |i), where i = ii, . . . , i„ and U denotes 
an arbitrary phase gate. We call a phase gate, U, pure 
phase gate of order k if it can be written as C/ = 1+ (e"^ — 
1) |1) (1| '. A phase gate will be said to be of order k 
if it can be decomposed into pure phase gates which arc 
maximally of order k. Note that any phase gate can be 
decomposed into pure phase gates. We will say that a 
qubit is overlapping if there exist at least two pure phase 
gates in the decomposition of [/, which act non-trivially 
on this qubit. If there are two or more qubits, say qubit 
1 and 2, for which there exists more than one pure phase 
gate, which act on them simultaneously, we say that the 



interaction possesses an overlapping edge between qubit 
1 and 2. 

Let us now consider an arbitrary bipartite splitting, 
A, B. where A {B) denotes a subsystem composed out 
of A^^ {Nb) qubits respectively. We write the n~qubit 
LMES, I*) as I*) = C/|+)®", where n = Na + Nb 
and U denotes a phase gate. Since U can be de- 
composed into pure commuting phase gates, we have 
U = UaUabUb, where Ua {Ub) is only acting on 
the qubits within A (B) respectively and Uab is act- 
ing non-trivially on both, A and B. Then, the re- 
duced state of system A can be written as pA ~ 
UAtTB{UAB{\+){+\)'^'-''^U\B)U\. Thus, we will focus 
in the following on the computation of the reduced state 
of I*) = Uab 1+)^" = Eie'"' |i). Writing |i) ^ |i^) jis) 



i{c 



<) 



De- 



wc have pA = Ei^j^,kB e"^"'A.«i5 -■>a.-b' \ia) {3a\ 
pending on the correlations of ai the reduced state can 
be computed efficiently, or not. We are going to con- 
sider now some simple examples of LMESs for which 
the reduced state can be computed efficiently. In gen- 
eral, however, this will not be possible, as will be ar- 
gued below. If 1^) is a product state, i.e. ai = J^k'^ik^ 
then PA = Y^iAJA e'^"'-^""-'-*^ ju) (Ja| is a pure product 
state. The second class in the hierarchy of LMESs (see 
Eq (2)) constitute the so-called weighted graph states 
(see e.g. [3, 4, 7]), which include (up to local uni- 
taries) the well-known stabilizer states [8]. There, ai 
is a polynomial of degree two in i, i.e. ai = i-^Fi 
(up to local unitaries), where the n x n matrix F = 
|0) (0| (8 Fa + |1) (1| ® Fb + (|0) (1| + |1) (0|) ® Tab is 
called adjacency matrix. In this case, the reduced state 
has the form (up to normalization) 



PA 



E e'^'^^^^'^'^^^^^^^^^^o,,^,,JsA){^A\, (3) 



where c?. =» 



n^1(l + e*^^^^''''^'-*"'"'^^^)> with e, de- 
noting the standard basis. Note that cj^.j' and therefore 
PA can be computed efficiently [9] . 

In general, however, it will not be possible to de- 
termined the reduced state of an arbitrary LMES effi- 
ciently. The reason for this is the following. Let us con- 
sider the n-qubit LMES, |*) = U\+)^"', where U is a. 
product of pure phase gates of order three, i.e. ai is 
a monomial of degree three in i. Let us now compute 
the reduced state of qubit 1. Without loss of general- 
ity we write U = ViVq, with Vq acting as the identity 



on the first qubit, and Vi = |0) (0|^ (g) 1 



11 



C/i, 



where Ui is a product of two-qubit phase gate. Thus, 



ID, c/i 



i{n-l) 



)■ 
n-1) 



we have \^) = Vo{\0),\+)"'^''-''> + \^,,u. 

Writing the weighted graph state, |0) = Ui 

as 10) = Es-e'*"^^*"!*) and noting that (""i)® (-h| ^ = 
^^TTTT Vs, we see that we would need to efficiently com- 



72" 



pute J2s^^'''^^ in order to efficiently compute pi. It is 
however very unlikely that this expression can be com- 
puted efficiently for an arbitrary F, since it resembles 



the partition function, which is in general not efSciently 
computable. 

We will use now the discussion above to understand for 
which states it will be possible to compute the reduced 
state efficiently. Suppose that F is block-diagonal, i.e. 
r = ^Fi, where the dimension of each of the block- 
diagonal matrices, F^ is independent of n. Then, it is 

easy to see that J2s^^^ ^"^ '^^^ ^^ computed efhciently 
and therefore the reduced state can be computed effi- 
ciently. Note that for the state |^) this means that qubit 
1 is independently coupled to several other qubits, but 
among the set of coupled qubits there is no coupling. 
That is, the number of overlapping qubits which are cou- 
pled to qubit 1 is small. Similarly, one can see which 
m~partite reduced states of LMES with higher order can 
be computed efficiently. 

Let us now compute the reduced state of certain 



:j3i32 



LMESs, 1^) 



t^l+> 



(gin 



which will become rele- 



vant in the subsequent sections. For a single qubit, 
which we choose here without loss of generality to 
be qubit 1, and an arbitrary interaction, U, which 
we write as before as U = ViVb, with Vq acting as 
the identity on the first qubit, and Vi = |0) (0|]^ ® 
1 + |1) (l|i ® Ui, where Ui = Ei e*"' |i> (i| we find 
pi = i(l + |0) (1| ("-!)« (+1 Ul |+)®("-i) + h.c), where 
(n-i)<» ^_^| ul |_^^)«^(n-l) ^ ^. e~*"' . Supposc now that 
Ul is a product of pure phase gates, Ui, each acting on 
ki — I different qubits with phases 4>i . Then, the reduced 
state is given by 

pi = 1(1 + |0) (1| n ^(2''-' - 1 + e-^'O + h.c.).ii) 

As another example we compute the single qubit re- 
duced state of the state U{(j)) |-|-) , where U{(f>) is a pure 
phase gate, i.e. [/(</>) = 1+ (e*"^ - 1) |1) (1|®". It is easy 
to show that pi{cj)) = tr^u buti('^(^) 1+) (+1" U{(j}Y) has 
the following form 



PiW 



2n-i 2""^-I + e 



-itp 



(5) 



Similarly the reduced state of several qubits can be 
determined. 

Another way to compute the reduced state is to 
make use of the notion of projected-pair-entangled state 
(PEPS) [10, 11]. The idea is to rewrite the quan- 
tum state by substituting physical qubits on which sev- 
eral phase gates are acting on by a number of virtual 
qubits. The transfer from the virtual system to the 
physical one is done by projectors acting on those vir- 
tual qubits. For example, if qubit i is part of k pure 
phase gates of arbitrary order, we would introduce k vir- 
tual qubits, one per phase gate. Then ik denotes the 
fcth virtual qubit of the physical qubit i. The following 
simple example of two three-qubit phase gates illustrates 



the idea: C/i23f^345 |+) 



^3'''t/l23iC/3.45 



\«i6 



P3^^^^= 72(10)3(0013^ 



.32 + 11)3 (ll|3i32)- The projector, 
P, is independent of the order of the phase gates applied 
to the physical qubit and independent of the phases. It 
only depends on the number of pure phase gates which 
are acting non-trivially on system 3. In general, the 
projector associated to a physical qubit i on which k 
pure phase gates are acting on has the following form 

Pr-'" = V2'"'(|0), (00 . . . 0|,^,...,^ + 11), (11 . . . 1|,^,...,J 
[19]. 

Let us now use the procedure described above to 
present the reduced state pA of an arbitrary LMES 



I*; 



AB 



= Uab\ 



\^{Na + Nb) 



Since the phase gates 



commute, the reduced state can be determined from 
an interaction pattern without the need of considering 
any time ordering. For each qubit in A and B, which 
interacts via several pure phase gates with other qubits, 
we introduce projectors and virtual qubits as explained 
above. Note that any pure phase gate is acting on 
different virtual qubits. That is, there is no virtual 
qubit on which two or more phase gates are acting 
on. All information about the coupling is then cap- 
tured by the projectors. An easy example for 2-qubit 
phase gates would be: C/i2C^23 1+) (+| ^^12^^23 = 

Pt''{Ui2,U2,3 1+) (+!' ul^ul,)Pt'^^ 



:)2i,22 



(Pl2i 



\pt2l,22 
P22 3J-f2 



where P, 



11)2(1112122) andp.j =C/y 



'U^. 
V 



V2{\0), (00|2^2, 

with (i,j) e {(l,2i), (22,3)}. In general, we 

define for any qubit bj in P, where '>Ti{bj) 

pure phase gates are acting on, the projector 



|l)fc,(ll 



1 



v^""(|o).(oo-;Ok,.....<.,+ 

) and use a similar definition 



with 



for any system Sj in A. Then the reduced state of 
system A can be written as pA — tr^d^) (^|) = 

PAtrB(PB <»i Pi (8){fc/} Pkl ■ ■ ■ <^{o....n} Po,...nPl)PA: whcrC 

Pb — ®j Pfc . ' and Pa is defined analogously. 

The operators p denote, like in the example above, pure 
LMESs. For instance, pi denotes a pure state which is 
obtained by applying a 1-qubit phase gates on the state 
1+) state, the operators pki result from 2-qubit phase 
gates apphed to |-l— 1-) and so forth. We will use this form 
later to compute the reduced state of special cases. 



III. DECOHERENCE MAPS 

We consider a system composed out of Ns system 
qubits, which will be denoted by si^ . . . sjsi^, which are 
interacting with Nb bath qubits, which will be denoted 
by 61 , . . . biVg . The qubits interact via some arbitrary 
phase gate, which will be denoted by Usb- The evolu- 
tion of the system qubits is governed by the completely 

positive map, £{ps) = tTB{UsBPs<^\+) (-t-l**^^ C^ss)- 1^^ 
this section we will first of all consider the most general 
dephasing map and show how it can be represented in 
terms of the reduced state of some LMES. After that we 



will represent the dephasing map in terms of Pauli opera- 
tors. Since we are interested in the decoherence caused by 
such an interaction we will restrict ourselves in Sec. Ill A 
to those interactions where non of the system qubits is 
directly interacting with another system qubit. As men- 
tioned before, we will call those maps purely dephasing 
maps. We will show that all of them are separable and 
derive a simple expression for them. Then we will inves- 
tigate how the correlations among the bath qubits will 
affect the evolution of the system qubits. 

As before, we decompose the general phase gate, UsB^ 
into UsB = UbUsbUs, where Ub {Us) is acting ex- 
clusively on the bath (system) qubits respectively and 
UsB describes the interaction between the system and 
bath qubits. Thus, we have £{ps) = UstTB(UsBPs <8i 
|-|-) (-|-| ^ Ul.g)Ug. In order to understand how de- 
coherence can affect the properties of the system state 
it is therefore sufficient to consider the map £{p) = 
tTBiUsBPs (» 1+) (-h|®^^ uIj^). In the following we will 
focus on this expression. 

Note that Usb can always be written as Usb — 
Sk \^)b (^I ® ^k, where [/k denotes a phase gate acting 
only on the system qubits [20] . Using this decomposition 
it is easy to see that 



£{ps)^J2^kPsUl^asQ 



PS, 



(6) 



,tS)Ns J A 



U^ and denotes the 



where as = Z^k'^kK 
Hadamard product [12]. Note that as = f ([+) (-h[®^^) 
is the (unnormalized) reduced state of the LMES 
Usb 1+) ", with n = Ns + Nb- In order to obtain the 
second equality in Eq. (6) we made use of the fact that for 
any n-qubit state, p and any phase gate, U, the following 
holds UpW = [U \+} {+f W] O p. Maps which can be 
written as in Eq (6) are called diagonal or Hadamard 
maps. For these maps several interesting properties in 
the context of optimal output entropy and purity have 
been recently shown [13]. 

It is important to note that the evolution of the sys- 
tem is completely determined by the reduced state as- 
Thus, all the correlations which are established between 
the system and the bath during the evolution are com- 
pletely captured by the reduced state of the correspond- 
ing LMES, Usb |+) "• As mentioned above, the com- 
putation of this reduced state gets more and more com- 
plicated as the number of overlapping edges, or qubits 
increases. Physically this corresponds to the situation 
where more and more bath qubits take part in joint col- 
lisions with system qubits. 

Another way of representing the dephasing map is 
to write it in terms of Pauli operators. To do so, we 
recall the so-called Choi-Jamiolkowski isomorphism 
between completely positive maps and positive semidef- 
inite operators (states). For a completely positive 
map, £ : B{Hi) — > B{'Hi), the corresponding state, 
£■12, is acting on the Hilbert space Hi Hi- It is 
defined as E12 = (Si h){Pi2), where P = |$+) ($+| 



with |$+) = zf™^^'^|M). On the other hand, 
we have, £{p) oc tr23(i?i2/53-P23)- In the case of de- 
phasing maps, the corresponding Choi-Jamiolkowski 
state, E, has the simple form i?si,...sjvc,oi,.aNs = 



£s,,...s., ® la,....a„,({8);i^ll</>+),^,, (0+1) 

{8),A.a.(^s ^ \+){+f''')<S)^PU, where 

p^^^^ = v^(|00)(00| + |11)(11|) and a, denote the 
auxiliary systems. That is, the Choi-Jamiolkowski state 
E can be rewritten in terms of the reduced state as of 
the LMES, UsB \+f- 

We will now use the operator E to rewrite 
the map in terms of the Pauli operators. 
We decompose E in the Bell basis, E = 

En,...,.„.j,....,,„ Ki::::t l*^i, • ■ • , *^J (*.i, ■ • ■ , *.„!, 

where |*jj = a^^ (g) ]1|$+), with ik € {0, . . . , 3} and 
CTo = 1, cTi = ax, a2 = ay, as = az denote the Pauli 
operators. In order to present the map in a compact 
form we use the notation |^i) = [^i^, . . . , ^If^^). Using 
now that £{p) oc tr23(i?i2P3-P23) wc find 



£(p)«^AJ((g)a,Jp((g)a,J, 



(7) 



fc=i 



1=1 



with A; = (^i|i?|^j). Due to the projector P occur- 
ring in the expression of E, which projects onto the sub- 
space spanned by {|^o) , j'l's)} we find that A j = unless 
ik,ji G {0,3} for aU 1 <k,j <n. 

Let us now consider some simple examples. First of 
all we consider the simplest case where a single sys- 
tem qubit is coupled to several bath qubits via two- 
qubit phase gates. This scenario was already treated 
in [4]. There, the decoherence map for the system 
qubit, denoted here by 1, which interacts with some bath 
qubits bi G Bi is given by £{p) = trBdlhigSi Uibiipi 
\+)b,. (+1) Ub, Ulj = X°p + \\azpaz + Xlltpa, - a^pt), 
with A[S = (1 + rcos7)/2, Aj = (1 - rcos7)/2 and 
Ag = (irsin7)/2 where r = Hfc es cos(0i&;/2) and 
7 = X^b es('/'ifci/2)- The easiest example when consider- 
ing 3-qubit interactions describes the situation where a 
single 3-qubit phase gate is acting on one system particle 
and two bath qubits denoted by bi bj. We find 



£{p) = trB{Uib,b,ip® 
= XqP - Xlazpaz 



Ibibi 



wLb.) 



Xl{lpaz-azpt) 



(8) 



with Aq 



^(7 



i(l — cos((/))) and 



8V I cos(0)), Ai 
|sin(0). 

Investigating now a sequence of interactions with three- 
qubit phase gates forces us to consider different cases, 
namely the cases of independent collisions or collisions 
with overlapping qubits/edges. We will discuss the maps 
describing the different coupling scenarios in Sec. IV B. 
For two system qubits where each of them couples to 
the bath qubits via one three-qubit phase gate we end 
up with three different cases. When the bath qubits are 



non-overlapping wc have £{p) = £i® £2{p) with £i de- 
noting the single qubit map for a three qubit phase gate 
(see Eq. (8)). If the system qubit 1 is interacting with 
two bath qubits, say hi 62 and system qubit 2 with &2 and 
&3, i.e. 62 denotes an overlapping qubit we find E{p) = 

Aoo/o + Xllipdzt - (T,lp) + Xlliplcj, - ia,p) + 
\ll{(J,tpa,l) + \ll{la,pa,l + cr^lplcr^ - cr^cr^p - 



a^a^pla^) + \ll{a ^a ^ p(T zO z) with 

h 3cos((/)2) -f cos((/)i)(3 -I- COS(02))), 

- cos(02))sin((/)i), Xll = ^(3 -I- 

cos(0i))sin(02), AiSJ = ^(3 -H cos(<^2))sin2(0i/2), 



\ll{t(Ty_pCF^a:, 

•^00 = 32 (^"^ 

\30 

-^00 



32' 



\30 
-^03 



30 — 16 ^ 

^(sin(</)i)sin(02)), A^i = ^^(3 
cos(0i))sin2(<^2/2), \ll = j^(sin2((/,i/2)sin(02)), Agi = 
j^(sin(0i)sin2(02/2)), Ai = i(sin2(0i/2) sin2(02/2)). 
For the case where 63 coincides with 61, i.e. when there is 
an overlapping edge, the structure of the map is the same 
as in the overlapping qubit case. The coefficients however 
change to Aj^j] = ^(13 -|- cos((^2) + cos((/>i)(l -I- cos((/)2))) 



and 



\30 

^00 



r(cos^((/)2/2))sin(^i) and A 



O)sin(02) Xll 



icos2(02/2)sin2(0i/2) 



\ll = 3^(sin(</)i)sin(</.2)) A8i = i(cos2(0i/2)) sin2(</.2/2) 

Ksin(0i)sin2(</.2/2)) 



16 

Aig = |(sin2(0i/2)sin(<^2)) Ai 

Aii = i(sin2(,/,i/2)sin2(<^2/2)). 

In summary, wc have shown here three different ways 
of representing dephasing maps, namely in terms of the 
Kraus operators (see Eq. (6)), as a Hadamard product 
with the reduced state 175 (see Eq. (6)) and in terms 
of the Pauli operators (see Eq. (7)). In the following 
we will mainly use Eq. (6) to investigate the effects of 
decoherence. 



A. Purely dephasing maps 

Let us now focus on those interactions where each 
system qubit is only interacting with bath qubits, but 
not with another system qubit. We will call these maps 
purely dephasing maps, because, as we will show, they 
are all separable, i.e. they are of the from £{p) = 
Y,,P^A^ (S) B, (S) ■ ■ ■ (S) N^pA\ «) B| «) • • • O iVt. Note that 
those maps cannot generate entanglement and can be im- 
plemented (at least probabilistically) by local operations 
and classical communication. The reason, why we are 
especially interested in those maps is because we are in- 
terested in the decoherence effect of multipartite phase 
gates and not in generating entanglement among the sys- 
tem qubits. Apart from that, considering those maps 
allows us to make a fair comparison to the decoherence 
caused by two-qubit phase gates, which always leads to 
separable maps. 

In this section we will first show that any purely de- 
phasing map is separable. Then we will restrict ourselves 
to those situations where m bath qubits are interacting 
with all qubits in the system. The system qubits are 



also interacting with some additional, but not overlap- 
ping qubits. We will present a simple formula for those 
maps, which will allow us to study how bath correlations 
affect the evolution of the system. This kind of investi- 
gations will further allow us to determine the behavior of 
the system in the case of Markovian and non-Markovian 
interactions. 

As before we consider the map £{p) = tcsiUsBPs ® 
(1+) (+|)®^^t^5B)- In order to deal with the various in- 
teractions we are going to introduce the following nota- 
tion. Let / = {ii, . . . ,ifc}, denote a set of k indices for 
some k < Nb and Bj denote those k bath qubits, i.e. 
Bj = {&ij, . . . ,&i^.}. Then, Ub,,si denotes a pure phase 
gate acting on the bath qubits, Bj and on the system 
qubit si. Furthermore, I denotes the set of all sets of 
indices of bath qubits which are interacting with some 
system qubit, i.e. I = {I E {1, . . . ,Nb}^'', for some k < 
Nb s.t. 3; s.t. Ubi,si ^ !}• That is, for each I € I 
there exist |/| bath qubits Bj and the single system 
qubit, si which interact with each other via the non- 
trivial pure phase gate of order |/| + I, UB,,sr For 
each system qubit, I, Ti <Z T denotes the set of indices, 
/ = {ii, . . . , i/j}, for some k < Nb such that the bath 
qubits Bj ~ {bi-^ , . . . , 6^^. } interact with the system qubit 
I. That is, Ii = {I e {1,..., iVij}^'^, for some fc < 
Nb s.t. Ubj.si 7^ !}• Moreover, for each set of in- 
dices, /, we introduce the function d/(k) = Oie/^j- 
The pure phase gates Ub,,si can be written as Ub,.si = 
1-1- |I, . . . , I) (I, . . . , l\g ®{Ui^si —1) where C//,s, is a sin- 
gle qubit phase gate acting on the system qubit s/ [21]. 
Thus, we have Ub,,si = Ek 1^) (kl^^ ® uf'}^'' «) 1, where 



U 



di{k) _ 



jdj{k) _ 



j^^^ = 1 if di{k) = and c/;;;;^ = t//,,, if d/(k) = i. 



Since the global unitary, Usb can be written as a prod- 
uct of the unitaries Ubi,s,, i-e. Usb = 11/ D/ei, Ub,,s, 
we find 



Usb = E \^) (k| 



®t/^. 



(9) 



where 



U^ 



n ut 



d,{k) 



C/f (E) C/2 



'Uhs, (10) 



/ leXi 



with U^ = ri/ei ^/ s ■ "^^^ corresponding completely 
positive (CP) map is given by 



£iP)=j2u^p(^sy-^sQ 



p, 



(11) 



where all Ug given in Eq. (10) are local unitary phase 

gates and as = Ek t^s(l+) (+l)®^''(t^s)^ is a convex 
combination of product LMESs and is, in particular, sep- 
arable (see also Eq (14) below). It should be noted here 
that since the bath qubits can be in an entangled state, 
there is no reason, in general, why the map should be 
separable. However, in the case investigated here, all the 



interactions commute, which implies that instead of con- 
sidering an entangled initial state of the bath qubits, we 
could equally well consider the product state, 1+) ^ as 
initial state. This fact alone is however, not enough to 
conclude that the maps arc always separable, since the 
baths qubits couple to several system qubits. 

So far we have seen that, whenever each system qubit 
interacts only with bath qubits any decoherence map, £ 
can be expressed as a convex combination of local deco- 
herence maps, i.e. E(p) = J2iPi^iip)- Here, £i denotes a 
local unitary map on the system qubits. The structure 
of £i is determined by the number of qubits the phase 
gates are acting on and by the number of overlapping 
qubits/edges. As a next step we will investigate the in- 
fluence of the coupling between the different baths (for 
each system qubit) on the evolution of the system qubits. 
The more bath qubits interact with all system qubits the 
stronger the correlation between the baths (for each sys- 
tem qubit) is. We will derive for any of these correlations 
a simple expression for the corresponding map, which will 
allow us then to analyze this effect quantitatively. 

We denote by £m{p) the state of the system qubits 
after each of them interacted via a phase gate (not nec- 
essarily pure) of arbitrary order, n; -I- 1, with the bath 
qubits {&ij , . . . , bi^ } and some other ni — m bath qubits, 
which do not interact with any other system qubit. That 
is, £m denotes the map which results from an interaction 
UsB = Y\i ri/gi Ubi.si, where X is such that the follow- 
ing condition is fulfilled. There exists some fixed set of 
m indices, {ii, . . . , im}, such that 



V/el, {ii, 



i} C / and 



b) V / G I;, /' € Iv with I ^ I' it holds that / n /' = 

For instance, £o{p) will describe the system qubits af- 
ter all of them interacted with independent bath qubits. 
£i (p) describes a situation where all the gates acting on 
the system qubits are coupled to one bath qubit, etc.. 
Considering these kind of maps allows us to study how 
the effect of decoherence changes with the degree of de- 
pendency of the environments of the system qubits. For 
instance, how does the decoherence affect the behav- 
ior of the system if all system qubits interact with the 
same bath, compared to the one, where each system in- 
teracts with its own bath. In order to study all those 
cases we arc going to derive now a simple expression for 
the corresponding maps. To this end we introduce the 
following notation. Ji ^ {J & {1,...,Nb}'' s.t. 3/ e 
Xis.t.J U {I, . . . ,in} ~ I}, i.e. each J € Ji denotes 
the bath qubits, Bj which interact with the bath qubits 
{6i, . . . 6m} and the l-ih system qubit via a pure phase 
gate. J denotes the union of all sets Ji. Let J = 
{jij • ■ • JlJl} be some set of indices. For some vector k = 
[ki , . . . , kNB ) we denote by kj the | J|-dimensional vector 
with entries kj for j e J, i.e. kj ~ {kj^ , . . . kj^j^ ). More- 
over kj, will denote the set of entries of k, {kj^ , ■ ■ ■ , kj^,}, 
where for each jt with t E {1, f} there exists some J E J^i 



such that it G J. Using this notation we can now state 
the following theorem. 

Theorem 1. Let Usb = 11/ Il/gi ^Bi.s,, where Bi 
denotes \I\ bath qubits and si denotes a single system 
qubit. Moreover, let X be such that the conditions (a) 
and (b) are fulfilled for some set of indices, {ii, . . . , im}. 
Then the corresponding CP map, £ = tr^ [t/s_B(p 
|+)(+|«^«)C/t^] IS given by 



£{P) 



1 



2Nb 



((2 



Nb 



■)NB—m 



)p + £i(g)...(g,£Ns{p))l^2) 



where £i (a) 



Ek, uf,:"a{U^,:"V. Here, U^'^ = 



d./(k) 



Ujej^ UZj 



is a local unitary operation. 



Note that this implies that the number of bath qubits 
which are interacting with all system qubits, m, only 
changes the factor in front of p and the weight with which 
C/]^ is applied but leaves the rest unchanged. 

Proof. We assume without loss of generality that the 
set of bath qubits with which any system qubit is in- 
teracting, {&ij, . . . , bi^}, is {6i, . . . , bm} (condition (a)). 
Then, due to the fact that {!,..., m} C / V/ G 
I we have that d/(k) = D" i ^» OsG/Xii,...,™} ^» = 
c(k)nie/\{i,...,m}fc*. where c(k) = UZi^i is indepen- 
dent of /. Thus, we find for the local unitary trans- 
formations, Uj" in Eq (10), U^ = {Uj'y''^\ where 

}(k) 



u^ 



n 



/ex, 



Tjdrsii. 



Eq (11) and the fact that X]k:c(k)=o 



Using this expression in 

^ „ 2^B _ 2^3— m 

we find £{p) = 5^[(2^« - 2^B-™)p ^ £(^)]^ ^^ere 

£ip) = Ek=(.„+,....,fe„,)(8>/^^(8>/(^^)^- As men- 
tioned before, J^ denotes the set of indices of baths qubits 
which are interacting with the system qubits and the 



bath qubits 6i 



, bm- Hence, for each pair J E Ji and 



J' e Jv with I y^l' we have Jf| J' = 0- Thus, the sum in 
the expression of £{p) can be decomposed as a sum over 
all bath qubits with are interacting with the first system 
qubit, J7i, the second, J^2, and so on. Since non of the in- 
dices will occur twice, we have £{p) = fi (g) . . . (g) £ns (p) ) i 
with £i as defined above. 

D 

It should be noted here that this result is independent 
of the coupling between the bath qubits and a single sys- 
tem qubit. The influence of their couplings and orders is 
reflected in the details of the local maps £[, but does not 
change the general form of the map. 

In summary, wc have seen that whenever some (to) of 
the bath qubits are interacting with all system qubits 
via some arbitrary pure phase gates and if no other 
overlapping bath qubit exists, then the map corresponds 
to a convex combination of the identity (with a weight 
that depends on to) and a local map. The local map 
corresponds to the map one would get by condition- 
ing on the fact that all qubits bi, . . .bm are in the 



state 

(l+X+l)^ 

where 



|1). That is, £i{a) ^ trs(C/B,^^, 



,}-sU 






Ub, 



,},si' 



-bi,...6„ 



(l.---,l|n/ei,C^s/^^|l'--"'l>L...fc™- 
Let us now consider some examples. We wih assume 

that any system qubit, Z, is interacting with the bath 
qubits only via a single pure phase gate of order ri; + 1 and 
phase 4>i. Suppose further that the conditions (a) and (b) 
are fulfilled for some set of m indices, {ii, . . . ,im}- Let 
£m ''"' ^ denote the corresponding map. It is straight- 
forward to see that the maps £i in Eq (12) arc then given 
by ii{a) = 2"'-™[(l - 2™-"0<7 + 2--"'f/,((/)i)aC/i(</)z)t]. 
Here and in the following Ui{(j)i) denotes the single qubit 

phase gate with phase 0;. Using then that Y\i 2"'"™ = 
2NB-m ^g ^^^ 



-,ni,...,nArc 



ip) 



1 



Ns 



[{2"^~l)p + 6<)Ei{p)l (13) 



1=1 



where 8i{a) = (1 - 2'»-"0ct + 2"-"'C/,(0OCTt/z(0O^- 
For instance, if all system qubits interact with the 
same m bath qubits (extreme non-Markovian case, 
see Sec. IV A) we obtain £:;^-'"(p) = ^[(2^« - 
2^«-")p + 2^«-™(£:i ® ...® £ns){p)1 where £i[u) = 
Ui{(f>i)aUi{(j)iy). We are going to show in Sec. IV B how 
the effect of decoherence depends on the number of bath 
qubits which are interacting with all system qubits. That 
is, we will analyze there how the degree of correlations 
between the environment of the various system qubits is 
affecting the evolution of the system. 

Note that any map which results from a pure phase 
gate of arbitrary order can be reinterpreted as a mixture 
of the identity map and a map resulting from a two- 
qubit phase gate acting on the system and the bath qubit, 
which is prepared in the state |^) ~ ^/a |0)+Vl — a^ |1), 
for some appropriately chosen a > 0. That is, £i{(t) ~ 
{l-2^-^^)a+2"'-^'Ui{cj)i)aUi{(Piy = trsiC/ssd*) (^|® 
cr)f/g^], where Usb denotes the two-qubit pure phase 
gate with phase (j)i and a = 1 — 2™""' . 

Let us now also express the maps studied above in 
terms of the Hadamard product, £{p) = as Q p- First of 
all, we present the state as , for the situation when there 
is no interaction between the system qubits (see Eq (11)). 
We start with the LMES IV-) = Usb |+)^^+^«. As the 
next step we use the PEPS picture and introduce vir- 
tual qubits 6ij , . . . , bi^ for every bath qubit bi on which k 
pure phase gates are acting on. For the qubits which are 
non-overlapping, which means just one pure phase gate 

is acting on them, the projector P^^'^' "' '*" will simply be 
the identity and for the overlapping bath qubits it will 



have the form P^ , 



'K 



= V2' 



fc-i. 



|l)ft^ (11 . . . l|fc^ (,^ ). For / e Ii, BJ denotes the set of 
virtual bath qubits which interact with system qubit s/ 
via a phase gate Us, b^ ■ The LMES in the PEPS picture 

can be written as |V') = ® A'^' ' ''' 0; 0/ei I'fpsi.Bj) 
with IV'.sj.B^) = Us,B^ 1+) • As a first step wc 



trace over all qubits which arc non-overlapping. 

^bi,,....bi. 



We 



denote hy F ~ {i : P^'^' ' "'' = 1;,.} the set of non- 
overlapping qubits. Then we have as ~ trsdi/') (V'l) = 

trB\F[0,^Fn''"""''"' ®i P^„Br<S)^^Fn^f'"""'''1 with 

Psi,s- = trigF((8)/G/i IV'si.Bj) (V's,,By|)- Performing the 
trace over all coupled qubits leads to a summation over 
all operators ®/ P^j where p^^ = (k|ps,,Bv |k), with |k) 
denoting the computational basis of the overlapping bath 
qubits. We end up with 



as 



E^ 



P% 



(14) 



Examples of the reduced state in this form are given in 
Sec. IV. For the special case where all system qubits 
couple to the same m bath qubits via one pure phase gate 
(sec Eq. (13)) wc find ag = 4r[(2'" - 1)(|+) (+|)^^« + 

^^=W 2'"""0 1+) (+1 + 2™-"'C/, 1+) (+1 C//)][22]. As 
mentioned before, for any n-qubit phase gate, U it holds 
that U \+) i+f" WQp = UpW. Thus, writing as as 
as ~ X^k ^k 1+) (+1*^" C/j[ as we did here, leads directly 
to the map £(p) = X^k UkpU^. Hence, £ can be directly 
read off given as and visa versa. Depending on which 
property of the output state one is interested in, one, or 
the other way of looking at the map is better suited. 



IV. COMPARISON OF DIFFERENT 
INTERACTIONS 

In this section we will use the result derived above in 
order to compare different decoherence models with each 
other. Due to the large variety of different interactions 
there is of course no way to make a complete comparison. 
However, we will focus here on certain relevant aspects 
of the dependency of the evolution of the system qubits 
on various couplings. In Sec IV A and Sec IV B we in- 
vestigate how the decoherence depends on the degree of 
bath-correlations. In Sec IV A, we will compare the case 
where there is no dependency, that is, each qubit inter- 
acts with an independent bath, and all pure phase gates 
are applied to fresh bath qubits and the system qubits 
(Markovian case) to the other extreme case, where the 
system qubits are always interacting with the same bath 
qubits (extreme non-Markovian case). We will deter- 
mine the coherence time for an individual system qubit 
and compare the two extreme cases for different orders 
of the applied phase gates. In Sec IV B we will first de- 
termine the evolution of a single system qubit and will 
investigate how the decoherence depends on the degree 
of correlation of the baths. In contrast to Sec IV A we 
keep the order of the phase gate fixed and vary the num- 
ber of qubits which are interacting more than once with 
the system qubit. Then we will consider 2-qubit states 
and certain n-qubit states and will show how the evo- 
lution of the entanglement shared by the system qubits 
depends on the coupling within the bath. Due to Eq. 



(12) one might expect that the system gets less disturbed 
the larger m is. However the examples investigated here 
will show that this is not the case. In fact we will show 
that for certain states the entanglement is more robust 
against decoherence induced by individual baths {m = 0) 
than against the one where m is very large. In Sec IV C 
we will discuss the effect of the order of the applied phase 
gates on the evolution of the system. 



A. Markovian scenario compared to (extreme) 
non— Markovian scenario 

In the following we will call an interaction, Usb, 
Markovian, if non of the qubits is interacting more than 
once with any other qubit via a pure phase gate, i.e. 
V/, I' E I, I n I' = 0. We will compare the decoherence 
induced by a Markovian process to the one induced by a 
non-Markovian process, where some bath qubits interact 
more than once with some system qubits. This compar- 
ison will be performed for several different orders of the 
applied phase gates. 

In the Markovian case we have that each system qubit 
interacts independently from the other system qubits 
with its own bath. Thus, we have that £{p) ^ £i ^ 
...(81 £ns{p)- Ii^ order to gain some insight into how 
this coupling affects the evolution of the system we study 
how the local maps, £i acts on a single qubit. The sin- 
gle qubit is, at each step, interacting with rii different 
bath qubits. We consider here the case where the sys- 
tem qubit is interacting k times via a (ni + l)-qubit 
phase gate with its environment. The phases of the 
phase gates are denoted by (pi, with i = 1, . . . , fc. The 
evolution is governed hj £{p) = £k ° ■ ■ ■ ° £i{p), where 
£t{p) = aip + PiUi{(t)i)pUi{(t)i)\ with [/, a single qubit 
phase gate with phase 4>i and Ui — 1 — 2^"^ , /3i ~ 2""^ . 
The relevant parameter to be considered here is the off- 
diagonal element of p, poi ~ (0| p |1), which allows us to 
determine the coherence time of the system [23]. Due 
to Eq. (6) it is easy to see that after the evolution, 
the off- diagonal element of the output state will be 
PoiCg^, where a^g — {Q\crs |1)- For the sake of simplic- 
ity we chose now (pi = (p ^i e {l,...,fc}. We obtain 



<^s = £{\+){+\) = ^i'(l+)(+l) = Oz=i'^^, where al 
corresponds to £i and is given in Eq. (5) with ni = n— 1. 
Thus, the absolute value of poi will be multiplied by the 



factor 



1 + 



Zl C0S((/)) 



2"l Y" ' 2"i-l "'""^'^^ ' (2"l-l)2 

Let us now investigate how this factor and therefore 
the coherence time depends on the order of the applied 
phase gates. The case where ni = 1, i.e. the inter- 
action is described by 2-qubit phase gates, has been 
studied in [4]. There, \cr^s'\ ~ [cos((?!)/2)]'^, which can 
be approximated by e"'"'"^ /^ for small values of (f^k has 
been obtained. In case the system qubit interacts via 
three-qubit phase gates with its bath we have \(j'^g\ ~ 
[ii/lO-l- 6cos(0)]'^ « e-rn^'t' , For four-qubit phase 
gates we get \ct%^\ = [i y'SO -h 14cos(0)]'^ « e-^^-^'/i^s. 



The coherence in the Markovian process for a single qubit 
falls therefore off exponentially with respect to the num- 
ber of occurred collisions, fc. The coherence time in- 
creases as the order of the phase gates get larger. In 
particular the ratio of the coherence times between the 
3-qubit phase gate case and the 2-qubit phase gate case 
is — g— « 1,3 and between the 4-qubit phase gate case 
and the 2-qubit phase gate case is — g^ « 2, 3. 

When the gates are coupled, i.e. the bath qubits inter- 
act more than once with the system qubit, we end up with 
a non-Markovian process, since the processes are not in- 
dependent of each other. In Sec IV B we will investigate 
how the evolution depends on the degree of correlation in 
the bath for a fixed order of applied phase gates. Here, 
we will consider the extreme case, where the system qubit 
is interacting always with the same ni bath qubits. It is 
easy to see that in this case the evolution is governed by 
the completely positive map £'(p) = pi{(pi + . . . + (pn) Q Pi 
where pi{(p) is given in Eq. (5). There, the coher- 
ence term, |poi|, gets multiplied by the factor \(7g^\ = 

^i^ Y 1 + 2^rr=T cos(X]Li ^i) + (2"i-i)^ - Like in the 
Markovian case we investigate how the coherence time 
depends on the order of the applied phase gate. Again, 
in order to make a simple comparison we choose (pi = (p 
Vi. The 2-qubit phase gate case was already discussed in 
[4] where \u%^\ = cos(k(p/2) « (,-k^<t>'^/^ has been found. 
The approximation is valid for small values of k(p. Since 
in the non-Markovian process the sum of all phases ap- 
pears in a cosine we get an oscillating behavior for large 
values of fc. For the three- and four-qubit interactions 

g-3fc^0V32 and 



we find \(j^g 



[iv/10 + 6cos(fc(/))] 



\af\ = [iv'50 + 14cos(fc</))] « e"^'' "^ /^^s^ respectively. 
The coherence time is longer if we go to collisions with 
phase gates of higher order. This is similar to the Marko- 
vian case. Within the range of the approximations the 
type of decay of the coherence in the non-Markovian case 
is Gaussian whereas the decay in the Markovian case is 
exponential with respect to the number of collisions fc. 
Note that the coefficients in the exponent coincide. 

In summary, we see that for both cases, the Markovian 
and the extremely non-Markovian case, the coherence 
time increases with the order of the applied phase gates. 
This can be understood by noting that the higher the or- 
der is the more qubits get entangled to the system qubit 
during the interaction. Since this entanglement is gen- 
erated by pure phase gates, the entanglement between a 
single qubit and the rest gets smaller. More precisely, if 
the single system qubit is interacting with its bath via 
a pure phase gate of order ni -I- 1, we have £„, (p) = pi, 
where pi is given in Eq. (5). It is easy to see that if rii = 1 
and (p = n, the system qubit is maximally entangled to 
its bath. However, for (ni + l)-qubit phase gates with 
ni > 1, this is no longer possible since the off-diagonal 
termof£(p), 2-("i+i)(-l-h2"i -He-'"^), is non-vanishing 
for any value of 0. As mentioned already, the only dif- 
ference between the two extreme cases investigated here 
is the dependency on the number of collisions (in the 
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region, where the approximations arc vahd). Whereas 
the coherence term in the Markovian case decays expo- 
nentially with the number of collisions it does so with a 
Gaussian function in the extreme non-Markovian case. 



B. Comparing different bath dependencies 



In this section we discuss how the induced decoherence 
depends on how strongly the baths for the individual sys- 
tem qubits are coupled with each other. To this end wc 
keep the order, rii 4- 1, of the applied pure phase gates 
fixed and the same for each system qubit. That is, each 
system qubit is interacting with rii bath qubits. Out of 
these ni baths qubits m of them interact with all sys- 
tem qubits. Using the notation from before, we have 
V/, /' e I / n /' = {6ij , . . . , bi^ } for some fixed set of m 
bath qubits, {bi-^, . . . , &i„, }. The aim of this subsection is 
to investigate how the evolution of the system changes as 
a function of m, the degree of bath-correlation. We will 
first work out in detail how the decoherence affects the 
evolution of a single and two system qubits. Then we will 
consider the general case of Ns system qubits and inves- 
tigate the influence of the degree of baths-correlation on 
the evolution. That is, we will use Eq. (12) and Eq. (13) 
to determine the evolution of the system for the various 
scenarios. For a single system qubits we consider the sit- 
uation where the system is interacting with the bath via 
several phase gates. The number of overlapping qubits, 
771, will denote in this case the number of qubits which are 
coupled to the system qubit in each phase gate. Physi- 
cally this means that those m bath qubits take part in 
any collision with the system qubit. We will investigate 
here how the decoherence of the system qubit depends on 
this coupling. For more than a single system qubit we in- 
vestigate the scenario where each system qubit interacts 
with its environment via a single pure phase gate. There, 
m denotes, as explained above, the number of bath qubits 
which interact with all system qubits. We will also inves- 
tigate the evolution of entanglement shared between the 
system qubits. 

As mentioned before, depending on which property of 
the evolved state one is interested in, the presentation of 
the decoherence map in terms of the Kraus operators, or 
in the form £{p) — crs Q p is better suited. This is why 
we will give both representations in the examples inves- 
tigated here. Whenever we talk about the correspond- 
ing reduced state in this context, we refer to as, which 
denotes the reduced state of the LMES corresponding to 
the evolution (see Sec. III). The phase gates, acting on n 



qubits will be denoted by C/"((^i) = 1— (1 - 



=#i 



)|1)(1 



(gin 



and Ui ~ U^{4>i) will denote the single qubit phase gates 
with phase (pi. 



1. A single system qubit 

We consider a single system qubit which is interact- 
ing with its environment via k different (ni + l)-qubit 
phase gate. We denote by m the number of bath qubits, 
{bi, . . . ,bm}, which take part in all interactions with 
the system qubit. That is, V/,/' G I with I ^ I' 
it holds that 10 1' = {bi, . . . ,bm}- Then we have 
£{p) = 2-"'[(2"' - l)p + o'^^^£,{p)], where o denotes the 
composition and £,{p) = 25;^((2"i-'" _ l)p+u^pU|). 
Let us now explicitly determine the evolution as a func- 
tion of m. In order to see the effect of the coupling within 
the bath it suffices to consider only two different phase 
gates. We start with the Markovian case, where the phase 
gates are acting on independent bath qubits. The corre- 
sponding decoherence map is £'(p) = as Q p = £2° £1 (p), 
where £j{p) = XqoP + ^mUjpUJ with Aqo = 1 - 2""^, 
Aqi = 2^"i. Note that the reduced state, (T5, is given by 
as = Pi Q P2 where pi is the single qubit reduced state 
of [/"i+H-^p 1+)^"^+' (see Sec II A). 
Next, we investigate the scenario in which one bath 
qubit is overlapping. The map can be written as 
£[p) = as Q p = ^[p + £2 ° £i{p)], where £j{p) = 
Aiop + XiiUjpU] with Aio = 1 - 2i-"i and An = 
2i-ni ^ rpjjg reduced state for this case is given by 
^s = 2EL0P1 pI with pO = i 1+) (+1 and p\ = 
_i_((2"i-i _ 1) 1+) (+1 + u, 1+) (+1 [//), for z e {1, 2}. 
The decoherence map for the case of two overlapping 
qubits is given by £{p) = as Q p = \[^p -^ £2 ° £i{p)], 



2-ni 



where £j{p) = A20P + A2iC/jpC/| and A20 = 1 — 2 

The corresponding reduced state is as = 



,00 



,10 



01 



A21 = 22-"i. 

4EL,.=oPi^'^0P^'^ withp^" = 5 |+; ;+| = p- = p 
andpii = ^{{2-^-^-l)\+){^ + UM){+\Ut),ior 
i € {1,2}. In the case where all qubits are overlapping, 
i.e. the gates act on the same bath qubits, we obtain 
£{p) = asQp = 2-"i[(2"i - l)p + UiU2pUlul\. The 
reduced state as is then the single qubit reduced state of 



?ni + l 



u 



111 + 1/ 



the state C/"i+i('^i)C/"'+^(02) |- 
4>2) 1+) 

It can be seen from the equation above that for a fixed 
order of the phase gates, ni, the coefficient in front of 
p gets larger as m increases. Thus, the system gets less 
disturbed if many of the bath qubits are overlapping. 
In order to investigate how the coherence for the single 
system qubit is affected by the number of coupled qubits, 
TO, we plot the absolute value of the off-diagonal element 
of as for different values of m. Fig IV BI shows the 
value of \a]f\ for 5-qubit phase gates for different number 
of coupled qubits m = 0,1,3,4 and for (j)! = (t>2 — 4>- 
The extreme cases where no bath qubits are overlapping, 
m = 0, and the case where all bath qubits are coupled, 
771 = 4, correspond to the Markovian and extreme non- 
Markovian case as already studied in Sec IV A. One main 
difference between those two cases is that the the non- 
Markovian case shows an oscillating behavior in contrast 
to the Markovian one. When both phases are tt, the value 
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a, 1.00 




FIG. 1: For a 5-qubit phase gate the absolute value of 
(0| crs |1) is plotted for different values of m. The continu- 
ous line shows the Markovian case where the two phase gates 
have no overlapping qubit, i.e. m = 0. The dashed-dotted 
line corredsponds to one overlapping qubit {m = 1) and the 
dashed line to three overlapping qubits (m = 3). The dot- 
ted line refers to the extreme non-Markovian case with four 
overlapping qubits (m = 4). 



of \(7g^\ is maximal for the non-Markovian case, since 
there the map acts as the identity map. However, it is 
easy to see that in the Markovian case the value of |crc;°| 
is minimal then. From the cases m = 1 and to = 3 one 
can infer the transition between the extreme cases. The 
m = 1 case looks very like the Markovian case, whereas 
the TO = 3 case shows already an oscillating behavior. For 
small phases the Markovian case is less damped than the 
non-Markovian one. This can be explained as follows. In 
the Markovian case the value of \(j]P\ is multiplied k times 
with itself, where k denotes the number of phase gates 
applied to the qubit. Since k equals only 2 in our case, 
the value remains relatively large since \ag^\ w 1. The 
value of \<Jg^\ for the non-Markovian case is decreasing, 
because we add the two phases of the two applied gates 
and therefore the value of the cosine appearing in the 
formula gets smaller. This effect is for small phases larger 
than in the Markovian scenario. 



2. Two system qubits 

In this subsection we consider two system qubits which 
are both interacting via a {ni + l)-qubit phase gate with 
the environment. Similarly to the single qubit case we in- 
vestigate here how the degree of dependency in the baths 
influences the evolution of the system. Here, we will also 
study how the entanglement between the system qubits 
is affected by the different couplings within the bath. 

The decohcrcncc map for the case where two inde- 
pendent (ni + l)-qubit phase gates are acting on a two 
qubit initial state (m = 0) has the following form: 
£{p) = as Q p = £2 ® £i{p) with £j{(t) = (1 - 
2^"i)cr + 2~"i[/jcrt/j). The reduced state is given by 
as = pi ® P2 with Pi being the single qubit reduced 



i.e. one bath qubit is overlapping, the decoherence map 
is given by £{p) = as Q p = \p + ^£2 ® £i{,p) with 
£j{p) = (1 - 2i-"i)cr + 2^-^WjaU]. The corresponding 
reduced state is 0-5 = 2 X]a;=o Pi®P2^ with p'^ = i |+) (+| 

and pi = 2ir((2"^-^ - 1) 1+) (+1 + U. 1+) (+1 U}). The 
case where two qubits are overlapping, i.e. to = 2, is 
described hy £{p) ^ as Q p ~ j[ip + £2 ® £i{p)] with 
£j{a) = (1 - 22-»i)cr + 22-"iC/jcr[/j). The correspond- 
ing reduced state is given by crs = 4 ^j. f. ^q p[^ "^ ® 



'^^'^^ with p™ ^ i 1+) (+1 = pi° = p°i and "" 



P2 



4 1+/ {+\ = P = P and pj - 
1) 1+) (+1 + C^j 1+) (+1 U]). In the extreme 
non-Markovian case where both system qubits interact 
with the same ni bath qubits, the map can be written as 



£(p) = (l-2-»i)p- 



tUt^ 



2-"^UiU2pUlU^ with the reduced 
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(see Sec HA). For to = 1, 



ilar to the single qubit case investigated before, we see 
that the factor in front of p increases as to increases. 
Thus, the system gets less distorted the more bath qubits 
are overlapping. 

Let us now study the evolution of the entanglement 
between the system qubits for several different coupling 
schemes. As examples we consider ni = 2, that is the 
system qubits are interacting with the environment via 
3-qubit phase gates. We compare three different types 
of coupling scenarios: 1) the case where each gate acts 
independently on one of the system qubits (to = 0); 2) 
the case where one bath qubits overlaps (to = 1 = ni — 1) 
and; 3) the case where two bath qubits are overlapping 
{m = 2 = ni). We compute the entanglement of forma- 
tion of the output state £{p) for the different couplings. 
The entanglement of formation is numerically optimized 
with respect to a maximally entangled input state. Fig. 
IV B 2 shows the entanglement of formation of the op- 
timized input state as a function of the two phases, (f)i 
and 02, of the three-qubit phase gates. In case 1), where 
771 = 0, the minimum at (/)i = 02 = "■, can be explained 
as follows. As explained in Sec IV A considering a single 
qubit interacting with its bath via a pure (ni -I- l)-qubit 
phase gate with phase generates a maximally entan- 
gled state only if = tt and ni = 1. For larger values 
of ni, the entanglement is still optimal, however it is no 
longer maximal, for cf) = tt. Thus, if both system qubits 
are interacting with their independent baths, the entan- 
glement between them is minimized if both phases are ir 
due to entanglement monogamy. 

In the second case, it can be seen in Fig IV B 2 that 
the entanglement is minimized if one of the phases is tt. 
This fact might be explained by the following observa- 
tion. The map which is applied to the input state is given 
by £{p) = \p+ \£i ® £2{p) with £j{a) = \a + ^UjaU] . 
As mentioned in Sec HI A we can interpret £j (a) as the 
map corresponding to the interaction of the system qubit 
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FIG. 2: entanglement of formation optimized with respect to 
a maximally entangled input state for three different coupling 
scenarios a) two three qubit phase gates C/ibi 62^^263 64 acting 
independently on the two system qubits 1 and 2 b) one bath 
qubit is overlapping within the phase gates 1/11,11,2^26263 c) 
the phase gates Uibib2U2bib2 have one overlapping edge 



with a single bath qubit via a two-qubit phase gate, [/, 
i.e. Sj{a) = trs(f/|+) (+| ® aW). Using this fact it is 
now easier to see why there is a minimum if one of the 
phases is tt. Similar to the reasoning above the two qubit 
phase gates have its maximal entangling capability be- 
tween the system and the bath for = tt and therefore 
the entanglement between the system qubits after apply- 
ing the map £1 £2 is minimized due to entanglement 
monogamy. 

In the third case, where tti = 711 = 2, we obtain a max- 
imally entangled state whenever (pi = (j)2 or 4>i + (j)2 = 27r. 
This can be easily explained by noting that the output 
state is of the form £{p) = ap + (3Ui U2pUl ® Ul, for 
some values of a and /3. Now, if 0i =02, which implies 
that Ul = U2 the singlet state, |\E'-) = l/\/2(|01) + |10)) 
is left invariant under this map. Therefore, the output 
state is maximally entangled. In the other case, where 



(pi + (j)2 = 27r, which implies that U2 = Ul, |$+) is left 
invariant. Thus in both cases the optimal entanglement 
of formation, optimized with respect to the maximally 
entangled input state, is maximal. 



3. n system qubits 

We consider the evolution of Ns system qubits, where 
each system qubit. s/ is interacting with ni bath qubits 
via a pure {ni + l)-phase gate with phase (pi. The de- 
gree of bath-correlation is quantified by the number of 
overlapping bath qubits, to, which we again denote by 
61,..., 6m- Those bath qubits are interacting with all 
system qubits and no other bath qubit is interacting with 
more than one system qubit (see conditions (a) and (b)). 
The evolution of the system is governed by the completely 
positive map given in Eq (12). Let us now explicitly write 
down the completely positive map for the three cases: 1) 
m = ni, 2) m = ni — 1, and 3) m = 0. 1) For m = ni 
all system qubits are coupled to the same m bath qubits 
via a (ni -I- l)-phase gate, i.e. the baths of all system 
qubits are maximally overlapping. In this case we have 
£{p) = ap+P^.^U^pUl with a = l-2-"i and /? = 2~"i. 
2) For TO = rii — 1 all system qubits are connected to the 
bath qubits {&i, . . . &„i-i} and one additional bath qubit, 
which is not connected to any other qubit. In this case 
we have £{p) = ap + P{®£i)[p) with a = 1 - 2-"i+i 



and /3 = 2-"i+i. Here, £i{ct) = aia + jSiUiaUj with 
a/ = /3; = i . 3) For to = each system qubit is interact- 
ing with its own bath. In this case the evolution if gov- 
erned by the completely positive map £{p) = {^£i){p), 
where £i{cf) ~ aia + PiUiaU} with ai = \ — 2^"' and 
Pi = 2~"'. As can be seen above (see also Eq (13)) the 
decoherence map is always of the form £{p) = ap+j3£{p), 
where £{p) denotes a local map. Since the coefficient a 
is given by 1 — 2~™, one might expect that the system 
gets less disturbed the larger to gets, that is the stronger 
the coupling among the baths is. However, as we will see 
below, this is not the case in general. For instance, the 
entanglement of a state can be decreased faster when the 
system qubits arc interacting with a maximally depen- 
dent bath than it is if each system qubit is interacting 
with its own bath. 

As example we consider the input state |^) = 
^ |0)®^^ + yr^ |1)®^^ with p > 0. We assume that 
all interactions with the environment are equivalent, i.e. 
(pi = (p, V/. It is straightforward to compute the output 
state for the different couplings mentioned above. Let us 
compare here the two extreme cases, ?Ti = and m ~ ni. 
In order to quantify the entanglement of the output state 
we consider here the bipartite splitting of one system 
qubit versus the remaining Ns — 1 qubits. Since the in- 
put state and the interaction is symmetric with respect 
to particle exchange, we consider, without loss generality 
the entanglement, -E, between qubit 1 and the rest, which 
we measure with the log-negativity [14]. It can be easily 
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three-qubit phase gates, ni = 2, with q 
E^ndep > Edep iff [\{^ + ^cos{Tr/Ns)) 



shown that for the case, where each system qubit is in- 
teracting with its own bath we find Eindep = log2(l + 
2|v^V'r^^(l - 2-"i + e''^2-"i)^s|) whereas in the 
extreme non-Mar kovian case we have E^ep = log2(l + 
2|^yr^Vl - 2-"i +e^^s02-"i|). The couphng of 
each system qubit to its individual bath destroys less en- 
tanglement than a correlated bath if Eindep > Edep- We 
are going to show next that such a scenario is possible for 
an arbitrary number of system qubits, Ns- Considering 

-^ we find 

Ns 

< 1. Since 
the left hand side of this inequality can be easily upper 

bounded by [j] 2 this shows that in general it is not 
true that the entanglement is less distorted if the degree 
of bath-correlation is larger. 



C. Comparison order of gates 

In this section we investigate the effect of the order ni 
of the phase gates on the evolution of the system. We will 
fix the number of overlapping bath qubits m. The influ- 
ence of ni can be seen by looking at the decoherence map 
presented in Eq. (13), i.e. C^'''""^ (p) = ^[(2^« - 

2^B-m-^p + 2^^-^^0l\£i{p)], where Siia) = ((1 - 
2^-^')a + 2"'-'^'Ui{cj)i)aUi{(j)i)'<). The order of the phase 
gates, ni , does not affect the coeflficient in front of p, but 
changes only the weights in the maps £1. The larger ni 
the larger is the coefficient in front of a. Therefore the 
initial state is less altered by phase gates of higher order. 
This coincides with the fact that the probability of ap- 
plying a phase gate on the system qubit gets smaller the 
larger the order of the phase gate is. For example, if the 
system interacts with the bath via a 5-qubit phase gate, 
all four bath qubits have to be in the state |1) in order to 
apply the phase gate to the system qubit; for a 2-qubit 
gate, however, just one bath qubit has to be in the state 



V. GENERATING LMESS 

In this section we will consider the complementary pro- 
cess to the one investigated so far. We will use the cou- 
pling of the system qubits to the bath in order to generate 
a state of interest. This can then be used to either pre- 
pare the state using dissipation, or a unitary evolution. 

We call a complete set of commuting unitary observ- 
ables, {Ui}^^Q generalized stabilizer if the eigenvalues of 
Ui are ±1 and both, 1 and —1 are 2"~^-fold degenerate 
[24]. Furthermore, we call a n partite state, IvP) a gen- 
eralized stabilizer state, if there exists a generalized sta- 
bilizer, {U,}f^o, such that a, |$) = |$) Vi iff |$) = |^). 
Examples of generalized stabilizer states are all LMESs. 
For an arbitrary n-qubit LMES, |*) = Ul \+)^'\ the 



U^Hl 



where X denotes here and in the following the Pauli op- 
erator (Tj.. It is easy to verify that [/& 1$) = 1$) Vfc iff 
1$) = I'ii). Depending only on the phases ai, which 
define the LME, |^), the generators of the generalized 
stabilizer can be quasi-local, i.e. act non trivially on 
a small set of (neighboring) qubits. In this case, the 
methods developed in [15] can be employed to derive a 
quasi-local dissipative process for which the unique sta- 
tionary state is |\I'). Apart form that, one can also easily 
construct a frustration free Hamiltonian for which the 
unique ground-state is | ^) [5] . 

In the following we show how generalized stabilizers 
can be prepared by a unitary evolution. Therefore, we 
use n auxiliary qubits which are all prepared in the state 
|-|-). The unitary operation is composed out of control 
operations, where the auxiliary systems act as control 
qubits. The following theorem shows that a generalized 
stabilizer state can always be prepared by applying this 
unitary to an arbitrary input state of the system. 

Theorem 2. // |^) is a n-quhit generalized stabilizer 
state then there exists a complete set 0/ commuting uni- 
tary observables, {Ui}"^Q (the generalized stabilizer) and 
a set {Vi}2^Q with {Vi, Ui} = such that for all n-qubit 
states 1$) there exists a n-qubit state \(j)) such that 



m 



= u^u^ . . 



■ u^m 



xSSn 



(15) 



The control operations V? are of the form Uf = UifJi, 
where (/, = 1 ® |0)^ (0| + U^ ® |1)„ (1| and f7, = 1 ® 

\+)aS+\+v^®\-):s~v 

Thus, in order to prepare a generalized stabilizer state, 
one only has to prepare the auxiliary systems in the prod- 
uct state |-|-) " and apply certain control gates. Depend- 
ing on the properties of the state | *) these gates might be 
quasi-local, i.e. act only on a few (neighboring) qubits. 
Note that this theorem can also be used to derive a pro- 
cess using either dissipation, or a completely positive map 
in order to prepare the state, |^) [15]. 

Proof. First of all, note that the set {Ui} docs not 
only define the state [^I'), but a whole set, the com- 
mon eigenbasis of the commuting observables. Each ele- 
ment of this basis is uniquely defined by its eigenvalues. 
We use the notation |* 



ii,. ..,•»„/ 



generalized stabilizer are of the form Ut = UKXk{U^- "i^ 



'ph/ 



for the basis elements 
with eigenvalues ((-l)*i, . . . (-1)*"), with ik € {0,1}. 
Let us denote by S*^ = span{\^ki....M=o,...MJkie{o.i} 
{SI = spar7,{|*fci^...,fci=i,...,fc„)fe,g{o,i}) the range of 1-hf/i 
(1 — Ui) respectively. That is any state which is an 
cigenstate of Ui with eigenvalue 1 (—1) is within Sq 
(SI) respectively and vice versa. Note that both sub- 
spaces have dimension 2"~^. We define Vi as the unitary 
which exchanges Sq with SI- To be more precise, Vi ~ 

Sfci,...fc„e{0.1} \^ku..;ki=0,..;kj (*fci,...,fci = l,...,fc„| + h.c. 

Now, we decompose |$) into a state which is in Sq and 
one which is in 51\ i.e. |$) = (l+L/„) |'I') + (]l-C/„) |$) = 
|*n) + !*«)• Applying C7„ to the state |$) |-|-)^ leads 
to u^ m [+),,, = 1$^) (|0) + |1))„„ + 1$!) (|0) Hl)),„. 
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Applying now C/„ leads to U^ |$) |+)^^ = |$o) |+),^ + 

Vn\^i)\-)a„^ where both, |$°) and 14 |^i^) are ele- 
ments of ^o ■ Continuing in this way, we find that all the 
terms which occur are elements oi Sq f] Sq f] . . . f] Sq = 
{|^)}. This is due to the fact that the operators Vj only 
change a state in SI to a state in Sq and vice versa, but if 
the resulting state is an eigenstate of Uk with eigenvalue 
1 (—1), it remains an eigenstate of Uk with the same 
eigenvalue. 

D 



Note that this is a new method to generate two copies 
of an entangled state, where both of them are generated 
in one step. That is, the states are not generated 
independently. One might compute the entanglement of 
formation using this procedure and compare it to twice 
the entanglement of formation of a single copy of |^). 



VI. CONCLUSION AND OUTLOOK 



Theorem 2 can also be used to show how an arbi- 
trary state can be constructed using a generalized sta- 
bilizer state as a recourse. This can be seen as fol- 
lows. Due to Theorem 2 we have that for any n- 
qubit state, |<I>) there exists a n-qubit state \(j)) such 



that I*) 



U^Ui 



■ u^m 



\^n 



Thus, we have 



= (*L u^ui . 



E 



X 



(ii) 



U'iUi 



\®n 



Note that 






\^n 



(16) 



1*1, •■ 



Here, H denotes the Hadamard transformation, and 



X 



(0) 



U„ and Xf ' = V,{1 - U,) = (1 + U,)V,. 



4n) 



.(»„) 



Since Uj is hermitian, we have that (^| X^ '-'■■■ Xn" = 
(^?i,....i„ I- Denoting by U^ the unitary which trans- 
forms the computational basis into the basis |'I'ii,...,i„), 
i.e. U^ = I]ii,...i„ l*ii,. ..,»„) (Ji,..., Jn| we have 



H 



Onrrt 



uu^) 



(17) 



Thus, an arbitrary state 1$) can be prepared by apply- 
ing the controlled operations (C/f)''' to the state |^) |(/)), 
where |0) is given in Eq. (17). For any LMES we 
have 1^) = t/* |0, . . . ,0). As before, one can define 
a phase-gate, t/^, which is diagonal in the computa- 
tion basis, such that, |^) ~ U^h\+) "? and, more gen- 
erally, f/^, = Upf^H®"-. In this case we would have 
10) = (t/p* )t 1$). Choosing |$) = |-H)^" we obtain 
|0) = 1^*), where |^*) denotes the complex conjugation 
of 1^) in the computational basis. Hence we have. 



I*) I** 
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(18) 



It should be noted here that for an arbitrary number 
of qubits, n, there exists a n~qubit LMES, JVP), such 
that 1^) and |^*) are neither LU-equivalent, nor can 
be transformed into each other by more general oper- 
ations, namely local operations and classical communi- 
cation [16]. For Graph states, into which any stabilizer 
state can be transformed by local unitary operations, and 
also for certain LMESs the phase gate is hermitian, which 
implies that |^) = |^*). Therefore, we find for those 
states 



I*) I*) 
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U" 



\ ®n 



(19) 



We investigated the effect of decoherence induced by 
many-body interactions. In our decoherence model we 
considered multi-qubit Ising interactions between the sys- 
tem and bath qubits. These interactions could result 
from multipartite collisions. The corresponding unitary 
was modeled as a multipartite phase gate. Due to the 
large variety of different interaction patterns we focused 
here on some relevant aspects of decoherence. We mainly 
considered purely dephasing maps, where each system 
qubit is coupled to another system qubit only via the 
environment, i.e. there is no collision between system 
qubits. We showed that all purely dephasing maps are 
separable. In order to analyze the effect of correlations 
within the bath on the evolution of the system we consid- 
ered those situations where m bath qubits are connected 
to all system qubits, which can also interact with other, 
non-overlapping bath qubits. We derived a simple ex- 
pression for those maps and showed how the effect of 
decoherence depends on the number of commonly shared 
bath qubits, m. This enabled us to compare all the cases, 
from the one where each system is coupled to its own en- 
vironment to the one where all system qubits are coupled 
to the same environment. Furthermore, we compared the 
Markovian to the Non-Markovian processes. We also 
showed how the evolution of the system is influenced by 
the order of applied phase gates. Note that these results 
can also be generalized to higher dimensional systems. 
The coupling to the environment can not only destroy 
the coherence and entanglement in the system, but can 
also be used to accomplish certain tasks [15]. We showed 
here how the bath qubits can be used to prepare a de- 
sired state. In particular, we derived a unitary opera- 
tion, acting on the bath and system qubits which can be 
employed to generate certain multipartite states. This 
investigation also lead to a protocol which prepares two 
copies of certain multipartite states in a single step. 

An other interesting question in this context is how 
the entanglement established between the system and 
bath qubits is related to the induced decoherence. For 
instance, the dephasing map £{p) — as Qp is a com- 
pletely dephasing map, i.e. £{p) = po, where po is a 
diagonal matrix whose diagonal is the same as the one 
of p, for any state p, iff as = 1. That is the map de- 
stroys all the coherence in the system iff the correspond- 
ing LMES is maximally entangled between the system 
and the bath qubits. A simple example would be the 



state 1*1 



SB 



UsB 



\®2n 



\^+)sB- Here, |$+) de- 
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notes a maximally entangled state shared between a bath 
and a system qubit. Obviously, for states, p, whose diag- 
onal elements are all equal, such a map would completely 
depolarize the state. Examples of such states are arbi- 
trary convex combination of product LMES (all in the 
computational basis). 

An other interesting topic is the characterization of de- 
cohcrcncc free subsets, i.e. the characterization of those 
states which are left unchanged under certain dephas- 
ing map. For instance, if we consider Ng system qubits, 
which are all coupled to the same bath composed of 
m qubits, i.e. ni ~ m + 1 for all I we have (see Eq. 



t 

Ns- 



(13)) £{p) = ap + pUi® ...® UnsPUI ® ...®U^ 
Thus, a state p is left invariant under this coupling iff 
C/i (8) . . . (g) UnsPUi (8i ■ . . <8i t/j^g = p- Choosing the phases 
such that all local unitaries coincide, we see that the state 
p will be invariant iff C/®^sp([/»JVs)t = p. states which 
fulfill the above condition for any unitray U are called 
super-singlets [17]. Since we consider here only dephas- 



ing maps, the state must only be invariant under local 
phase gates, which implies that the decoherence-free sub- 
set contains more states than the super-singlets. For 
instance, in the case of two qubits, the so-called Werner 
states, which are of the form p = {l-p)/4:l+p\'$^) ('f^|, 
where 1 > p > and |^^) denotes the singlet state, 
characterize all states invariant under U ® U , for any 
unitary U [18]. However, the decoherence free subspace 
of the map described above contains any state of the form 
p = a |00) (00| + /3 111) (11| + 7 [*+> (*+[ + <5 [*-) (*-[, 
for a,/3,7, (5 > 0. 
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as I*) 
V2{\Q) • 
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(0123)(|^+)3^J<^+) 
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Note that in order to make the corresponding completely 
positive (CP) map trace preserving we would simply need 
to consider \/2^^Usb- In the following we will however 
not explicitly write this normalization factor. 
Note that any phase gate can be decomposed in this form. 
Note that this corresponds to a map, which is a convex 
combination of the identity map and the one which re- 
sults from 2-qubit phase gates, as it should be. 
As explained in the introduction, the phase gates cor- 
respond to a evolution due to the Hamiltonian of the 
form HsB = X^i^^i'^z- The time dependence of the uni- 
tary UsB ~ e^'^^^* is hidden in the phases of the phase 
gates. This is why we speak of coherence times, even 
though we deal with time independent gates. 
Note that for stabilizer states, the stabilizers Ui would 
all be elements of the Pauli group. In general however, 
there is no reason to restrict to those stabilizers. 



